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ABSTRACT 

We present results on spectral variability of two Galactic black hole X-ray binaries, 
GRO J1655-40 and GX 339-4, in the hard state. We confirm a transition in behaviour 
of the photon index with luminosity, such that the well known decrease in X-ray pho- 
ton index with decreasing luminosity only continues down to Lboi ~ O.OIXe- Below 
this point the photon index increases again. For Comptonisation models, this implies 
that the ratio of the Compton luminosity to seed photon luminosity, £h/^s, changes 
with bolometric luminosity, consistent with a scenario where seed photons change from 
cyclo-synchrotron at the lowest luminosities to those from a truncated disc. Alterna- 
tively, the transition could mark the point below which the non-thermal jet starts 
to dominate, or where reprocessed photons replace the viscous ones in an outflowing 
corona model. 

Key words: accretion, accretion discs - black hole physics - X-rays: binaries. 



O 



- 1—1 

x 



1 INTRODUCTION 

Galactic black hole binaries (GBHs) are strong sources of 
X-ray radiation. This radiation can be decomposed into two 
main components: soft thermal radiation originating in an 
accretion disc, and hard power-law like emission, probably 
originating in the process of Comptonisation of the disc pho- 
tons by hot electrons located in the so-called corona (see re- 
views by Remillard & McClintock 2006; Done, Gierlihski & 
Kubota 2007; and references therein). The disk and corona 
contribute at different levels to the total spectrum. Conse- 
quently, a number of different GBH spectral shapes have 
been observed and classified into two main spectral states. 
In the soft state the hot disc dominates the spectrum up 
to a few keV and is accompanied by a weak power-law tail 
with a soft photon index, F > 2. In the hard state the disc is 
much cooler and the total X-ray band is dominated by the 
radiation of the corona with typically a hard photon index, 
T < 2. 

The physical origin of the transitions between the spec- 
tral states is not clear yet. Nevertheless, phenomenologically 
the evolution of the X-ray spectral shape within a given 
state and during the transition from one state to another is 
relatively well understood. The GBHs in outburst track the 
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same, continuous q-shaped pattern in the so-called hardness- 
intensity diagram (HID, e.g. Fender, Belloni & Gallo 2004; 
Dunn et al. 2008). From this diagram it follows that ini- 
tially the X-ray spectra of GBHs are dominated by the hard 
colour and their intensity is low; the sources are in a hard 
state and reside in the bottom/right corner of the HID. With 
time, the intensity increases but the spectra remain domi- 
nated by the hard colour; the sources move upwards in the 
diagram. At some point, further increase in intensity results 
in the decrease of hardness, and the sources move toward the 
top/left part of the HID. Eventually, a transition to a soft 
state takes place, resulting in populating the top/left section 
of the HID. Towards the decay of the outburst the intensity 
decreases and a transition back to a hard state is observed. 
The soft-to-hard state transition occurs typically at lower 
intensity than the hard-to-soft state transition at the begin- 
ning of the outburst (e.g. Gierlihski & Newton 2006), which 
leads to a hysteresis pattern in the HID. Theoretical expla- 
nation of the hysteresis include effects of different source 
of coronal irradiation in the hard and soft states (Meyer- 
Hofmeister, Liu & Meyer 2005, 2009), and in the most ex- 
treme cases a non-steady state behaviour when the mass 
accretion rate changes by several orders of magnitude over 
a period of week, or shorter (Gladstone, Done & Gierlihski 
2007). The X-ray spectral evolution in the HID correlates 
with the X-ray variability properties (Belloni et al. 2005) 
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Table 1. Log of the RXTE observations. (1) The name of the source. (2-4) The year, and the date of the first/last observation we studied. 
(5) The hydrogen column density used in the spectral fits. (6—8) The distance (the uncertainty on the distance from the literature), 
inclination and black hole mass of the systems. (9) The number of analysed hard state observations. References for the entries in columns 
(5-8) are as follows: [1] Done & Gierlirlski (2003), [2] Orosz & Bailyn (1997), [3] Remillard & McClintock (2006), [4] Hynes et al. (2004), 
[5] Zdziarski et al. (2004), [6] Munoz-Darias, Casares & Martinez-Pais (2008), [7] since the inclination to GX 339-4 is unknown we assume 
30°, [8] Kolehmainen & Done (2010). 
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Figure 1. Hard state lightcurves of (a) GRO J1655-40 and (b) 
GX 339-4. See e.g. Done et al. (2007) for complete lightcurves. 
Bolometric luminosities in Eddington units have been calculated 
based on the best fitting models of Sobolewska et al. (2011, GRO 
J1655-40) and Model 2 in this work (GX 339-4) extrapolated to 
the 0.01-1000 keV band. 



and with properties of the radio emission from a jet (Fender 
et al. 2004). 

In the HID an increase in intensity is associated with 
an increase in the mass accretion rate, while a decrease in 
hardness corresponds to softening of the X-ray photon in- 
dex. A number of scenarios were presented to explain the 
spectral evolution of GBHs, i.e. the spectral softening (hard- 
ening) with increasing (decreasing) mass accretion rate (e.g. 
reviews in Remillard & McClintock 2006; Done et al. 2007). 
One of the scenarios is that of the truncated accretion disc 
(e.g. Esin et al. 1997). In this model at low mass accretion 
rates the standard accretion disc in the hard state is trun- 
cated far away from a black hole and a hot inner flow is 
formed close to the black hole. In such geometry, few disc 
photons can be intercepted by the hot flow and the result- 
ing spectra have hard photon index, characteristic of the 
hard state. As the mass accretion rate increases, the inner 
radius of the disc moves closer to the black hole resulting 
in an increased cooling of the electrons. Consequently, the 
spectrum softens and the source evolves to a soft state. The 
soft state spectra seem to be well understood in terms of 
the multicolour disc black body emission from an untrun- 



cated (extending to the innermost stable circular orbit) disc 
accompanied by a weak Comptonised hard tail. Alternative 
hard state models not involving a truncated disc were pro- 
posed (see Done et al. 2007 for strengths and weaknesses of 
each model). These include e.g. a patchy outflowing corona 
above a disc extending to the last stable orbit (Beloborodov 
1999a; Malzac et al. 2001) and the jet origin of the hard 
X-rays (Markoff, Falcke & Fender 2001). 

However, several recent studies have shown that in the 
hard state at very low mass accretion rates, below a few per 
cent of the Eddington rate, the spectra of GBHs begin to 
soften again while the luminosity decreases (e.g. Ebisawa et 
al. 1994; Revnivtsev, Trudolyubov & Borozdin 2000; Corbel 
et al. 2004; Jonker et al. 2004; Kalemci et al. 2005; Wu & 
Gu 2008; Dunn et al. 2010). Interestingly, a similar spec- 
tral evolution with luminosity has been reported also in the 
case of AGN. Above a certain luminosity, the photon index 
softens with increasing luminosity (e.g. Porquet et al. 2004; 
Shemmer et al. 2006; Saez et al. 2008; Sobolewska & Pa- 
padakis 2009), but at low accretion rates the opposite trend 
is observed: the photon index hardens while the luminosity 
increases (Constantin et al. 2009; Gu & Cao 2009), just like 
in GBHs. 

As mentioned, it is generally believed that the hard 
power-law like X-rays are produced from Compton up- 
scatter of soft disc photons by energetic electrons in a hot 
corona located close to the black hole. The high energy cut- 
offs observed in the hard state spectra suggest that the pop- 
ulation of electrons is thermal (e.g. Gierlihski et al. 1999; 
Rodriguez, Corbel & Tomsick 2003; Miyakawa et al. 2008; 
Joinet, Kalemci & Senziani 2008; Motta, Belloni & Homan 
2009). Within this framework, the main parameter that de- 
termines the spectral shape of the intrinsic hard X-ray con- 
tinuum in accreting sources is the ratio of heating-to-cooling 
compactnesses, £h/£ s , where compactness is defined as a di- 
mensionless luminosity, 



where L is the luminosity of a spherical region of radius R, 
and ot is the Thomson cross section. 

In this paper we study how the compactness ratio, 
ih/is, evolves with GBH luminosity in the hard state. We 
consider the rise and decay parts of the outbursts of two con- 
firmed black hole binaries, when their luminosity is less than 
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Figure 2. (a/d) The hard state evolution of the X-ray colour with the luminosity in Eddington units in GRO J1655-40 and GX 339-4; 
filled symbols — rise, opened symbols - decay. A hysteresis pattern can be seen in both sources, (b/e) The X-ray photon index as a 
function of the luminosity in Eddington units (the best fitting photon indices are those reported by Sobolewska et al. 2009). (c/f) The 
correlation between the X-ray colour and the photon index suggesting that in GRO J1655-40 and GX 339-4 the colour variations are 
caused by the photon index variations in the hard spectral state at all luminosity levels. The solid line represents the best fitting linear 
function; in GX 339-4 we fit only the decay data (see text). The observations that give only upper limits on the photon index were 
neglected during the fit. 



~ 0.1-0.3 of their Eddington limit. We fitted their spectra 
with a disc black-body component (to account for their disc 
emission in soft X-rays) and the eqpair model of Coppi 
(1999) to determine lh/£ s - Since this parameter is mainly 
determined by the disc/corona geometry, our aim is to infer 
and constrain the evolution of the source geometry in the 
hard state, and investigate which of the current theoretical 
models are consistent with our results. 



2 DATA SELECTION 

In this work, we re-analyse all the archival RXTE observa- 
tions of GRO J1655-40 and GX 339-4, available up to 2007. 
Both sources are well studied systems containing a dynam- 
ically confirmed black hole primary (see references in Tab. 
1). The GRO J1655-40 selected data cover in great detail 
the period of its 2005 outburst. In addition, the mass and 
distance to GRO J 1655-40 are relatively well constrained 
(Tab. 1; but see Foellmi 2009, and Sec. 4.3), which allows 
for accurate conversion of the observed fluxes into the Ed- 
dington luminosity ratios. The GX 339-4 selected data ex- 
tend over a time-scale of several years. They cover reason- 
ably well its 2002/2003 and 2004/2005 outbursts (o2 and 



o3, respectively) and include observations taken in the pe- 
riod of 1996-1999 (ol). Even though the mass and distance 
constraints of GX 339-4 are not as tight as those of GRO 
J1655-40, we choose this system for our study due to a large 
number of hard state observations displayed in its multiple 
outbursts. Both sources follow a well known q-shaped pat- 
tern in the colour-intensity diagram (e.g. Dunn et al. 2010), 
which means that their X-ray spectral properties are typical 
of GBHs. Hence, we believe that the data we selected are 
representative of a GBH behaviour during the luminosity 
rise and decay phases of an outburst. 

We extracted archival Standard 2 spectra for detector 
2, top layer, in the 3-20 keV band and HEXTE spectra from 
both detectors in the 20-200 keV band. The details on data 
reduction are given in Sobolewska, Gierlinski & Siemigi- 
nowska (2009). Standard methods were used to subtract the 
background, create response matrices and deal with system- 
atic errors. We obtained one spectrum per pointed observa- 
tion. 

Of all the available RXTE observations, we chose to 
study the data when the source's high energy spectrum has 
photon index T < 2, as found from the PCA+HEXTE 3- 
200 keV model fits including thermal Comptonisation com- 
ponent (Sobolewska et al. 2009). Figure [T] shows the X-ray 
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luminosity in Eddington units plotted as a function of time, 
for the selected RXTE observations of the two sources. The 
rising/decaying phases of the outbursts can be clearly seen. 
In all cases, the source luminosity was lower than 10% of 
the Eddington limit for GRO J 1655-40 and 30% for GX 
339-4. Consequently, both systems during these observations 
should be mainly in their hard state. Luminosity estimates 
are based on the Model 2 that we describe in the next sec- 
tion. In total, we studied 94 and 224 hard state observations 
of GRO J1655-40 and GX 339-4, respectively. 

In Tab. 1 we list the time period covered by the selected 
data, the hydrogen column density used in the spectral fits, 
the physical properties (i.e. distance, inclination and black 
hole mass) we adapted in this work, and the number of 
the hard state observations that we studied. The distance 
to GRO J1655-40 and its mass are known quite accurately 
(but see Foellmi 2009). However, the estimates for GX 339-4 
suffer from significant uncertainties, so in this work we ad- 
dopted the canonical values (for GBHs) of 10 M and 10 
kpc for this source. 



3 SPECTRAL MODELING 

The process of Comptonisation of the soft disc photons in 
the hot corona is widely accepted as the mechanism respon- 
sible for the bulk of the power-law like continuum in the 
X-ray spectra of accreting objects. We considered two dif- 
ferent models to fit the 3-200 keV band spectra of our two 
sources. We included a multicolour disc blackbody (diskbb 
in Xspec) in both of them to account for the disc emis- 
sion. As for the power-law continuum, we used two different 
Comptonisation models: the thcomp thermal Comptonisa- 
tion routine of Zdziarski, Johnson & Magdziarz (1996), and 
the eqpair code of Coppi (1999). In both models we fixed 
the temperature of the soft disc photons at fcTbb = 0.4 keV 
because the RXTE low energy bandpass does not provide 
enough coverage to measure properly the temperature of the 
cool hard state discs. We modeled the spectra using Xspec 
ver. 11.3.2 (Arnaud 1996). 



3.1 Estimating the photon index 

Sobolewska et al. (2009) modeled the 3-200 keV 
PCA/HEXTE data of several GBHs. They used a model 
consisting of diskbb, thcomp, a Gaussian line profile to 
model iron Ka features, and a smeared edge (smedge) to 
account for iron Ka absorption. The complete model was de- 
scribed in Xspec as constant*wabs*smedge*(gaussian 
+ diskbb + thcomp), hereafter Model 1. The Nh for the 
absorption component wabs was kept fixed to the Galac- 
tic column density listed in Tab. 1. A constant component 
was added to account for differences between the PCA and 
HEXTE normalisations (the constant was fixed at 1 for 
PCA and was left free to vary for HEXTE). We used the 
Sobolewska et al. (2009) results to study the evolution of 
the photon index in GRO J1655-40 and GX 339-4. 



3.2 Estimating the heating-to-cooling 
compactness ratio 

Sobolewska, Siemiginowska & Gierlihski (2011) modeled the 
hard state GRO J1655-40 observations of the 2005 out- 
burst replacing the THCOMP model with eqpair, a hybrid 
thermal/non-thermal Comptonisation model (Coppi 1999, 
Gierlihski et al. 1999). The main parameter determining the 
spectral shape in the eqpair model is the ratio between the 
compactness of seed photons, i s , and hot electrons, in- This 
ratio, In/is, defines the spectral shape of the Comptonised 
continuum. It is a physical parameter (as opposed to the 
photon index in Model 1) and it depends mostly on the 
geometry of the accretion flow. Typically, the hard state 
spectra are characterised by lulls S> 1, soft state spec- 
tra by it, I i s <C 1 and very high/intermediate state spectra 
by lh/is ~ 1. The complete model was defined in Xspec 
as constant*wabs(diskbb+eqpair), hereafter Model 2, 
where CONSTANT and WABS were defined as in Model 1. 

As described in Sobolewska et al. (2011), 8 of 94 hard 
state observations showed residuals around 6-7 keV reminis- 
cent of an iron line, and so the LAOR relativistic line model 
(Laor 1991) was included to improve these fits (eqpair does 
not account for the iron line emission). With this addition, 
Model 2 fits to all 94 hard state datasets of GRO J 1655-40 
resulted in a null hypothesis probability higher than 0.05. 

In this paper we applied the same model to the selected 
hard state observations of GX 339-4. We found that 48 of 224 
datasets required an addition of the iron line, and again we 
added the LAOR component to Model 2. Following this, the 
model provided a good fit (with a null hypothesis probability 
higher than 0.05) to 220 of the 224 datesets. We therefore 
concluded that Model 2 is a good description of the hard 
state spectra of this source as well. 

We use the results of Sobolewska et al. (2011) for GRO 
J1655-40 and our results for GX 339-4 to study the hard 
state evolution of In I Is- 



3.3 Estimating the bolometric luminosity 

The spectra we studied are characterised by a photon index 
lower than ~2. We assume that they originate in a process 
of thermal Comptonisation, and so they peak in the EFe 
representation at energies of the order of the temperature 
of the thermal electron. Indeed, many of our spectra show 
a roll over at energies of the order of 100 keV. In order to 
estimate the total X-ray flux, and ensure that we do not 
neglect any energy radiated above 100 keV, we extrapolated 
the best fitting Model 2 spectra, and estimated the flux in 
the 0.01-1000 keV band, Fo.oi-iooo- Since the majority of 
energy in GBHs is radiated in the X-ray band, we consid- 
ered Fo. 01-1000 to be representative of the bolometric flux. 
Assuming isotropic emission, we defined bolometric lumi- 
nosity as L b oi = 47r_D 2 Fo. 01-1000, for a black hole located at 
a distance D. We study the spectral evolution of our sources 
as a function of bolometric luminosity in Eddington units, 
L bo i/L E , where L E = 1.3 x 10 38 M/M ergs s _1 for a black 
hole of mass M. 
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Figure 3. The hard state hcating-to-cooling compactness ratio of the corona, lh/&a as a function of the luminosity in Eddington units 
for (a) GRO J1655-40 and (b-d) GX 339-4. Filled symbols - rise; open symbols - decay. Solid lines indicate the broken power law fits 
to the data. The light gray data points were neglected during the fit (see Sec. 4.3). In GX 339-4 the best fitting break luminosity is 
^crit.GRO = (0.0061 ± 0.0002)L E , and in GX 339-4 the average of the fits to the three outbursts is L crit ,GX = (0.011 ± 0.002)L E . 



4 THE SPECTRAL EVOLUTION IN THE 
HARD STATE 

4.1 X-ray colour evolution 

We first considered the evolution of the X-ray colour with 
source accretion rate in order to study the spectral variabil- 
ity of the two sources in a model independent way, without 
any a priori assumption about the shape of the X-ray con- 
tinuum. Following Belloni et al. (2005), we defined the X-ray 
colour as the ratio of the observed count rates in the hard, 
i.e. 6.3-10.5 keV, and soft, i.e. 3.8-6.3 keV, X-ray energy 
bands. Our results are shown in Figs.[2^/d. Note that, tra- 
ditionally, this kind of plots show a count rate as a function 
of X-ray colour. However, we chose to plot the Eddington lu- 
minosity ratio for a more straightforward comparison of the 
plots in Figs[2ji/d with the following plots in the present 
work. 

The figures show clearly that the X-ray spectra of GRO 
J1655-40 and GX 339-4 vary with luminosity during the 
hard state. The X-ray colour follows a well defined pattern, 
which is similar in both sources. It first increases, i.e. the 
hard X-rays become more dominant, with increasing lumi- 
nosity until Lboi/iE reaches the value of ~ 1-3%, above 
which the colour decreases. This value obviously depends 
on the adopted distance and black hole mass estimates for 
each source, however the inferred spectral evolution with lu- 
minosity cannot be attributed to any specific model assump- 
tions because the colour was calculated using the observed 
count rates. We indicate the rising and decaying phases of 
the outbursts with different symbols and notice the hystere- 
sis pattern in both sources. 

4.2 Photon index evolution 

Although the X-ray colour is a model-independent measure 
of the spectral shape of the sources, it is not easy to interpret 
its evolution with source luminosity, as colour variations can 
be caused either by a variable normalisation of the disc and 
Comptonisation components, and/or by intrinsic variations 
of the hard band photon index. For that reason, Figs. rj}j/e 
shows a plot of the best fitting X-ray photon index from 
Model 1 as a function of source luminosity, for the same 
observations that are plotted in Figs. [2ji/d. The F values 



plotted in these panels correspond to the Model 1 best fitting 
results of Sobolewska et al. (2009). 

The best fitting V values for the lowest flux spectra 
have rather large errors associated with them. In fact, in 
some cases (specially in GX 339-4) , the best fit value reached 
the lowest boundary allowed during the fit, set to V = 1.4 
(Sobolewska et al. 2009). In these cases, the photon index 
estimates should be considered as upper limits. These effects 
complicate the determination of the photon index evolution 
with luminosity, at the lowest flux states of the sources. 

Overall, we believe that, the evolution of T with 
Lboi/LE as shown in Figs. [5p/e, shows a pattern which 
is similar to that of X-ray colour evolution (Figs. (2^,/d). 
As the source luminosity increases, F becomes harder, until 
Lboi/L-E reaches the level of a few per cent, above which 
the spectrum softens with increasing luminosity. This result 
suggests that, to a large extend, the X-ray colour evolution 
is caused by intrinsic photon index variations when GBHs 
are in their hard state. To strengthen this conclusion we 
plot the photon index as the function of the X-ray colour in 
Figs. HJi/f- It is apparent that the two quantities correlate 
and their dependence can be modeled with a linear function. 
During the fit we neglect the data points being upper limits 
(but we include them in Figs. [2J;/f), and we use the average 
of the upper and lower F errorbars (from Model 1) as the er- 
ror on T. We obtain Tgro = (-T60T0.04) xC+(2.90±0.03), 
Tgx = (-1.6 ±0.2) x C+ (2.9 ±0.2), where C stands for the 
X-ray colour. The fits result in relatively high \ 2 /d.o.f. of 
3.8 and 2.3 for GRO J1655-40 and GX 339-4, respectively. 
Nevertheless, it is clear that the model describes rather well 
the observed anti-correlation. In the case of GX 339-4 we fit 
only the data points of the outburst decay. The observations 
taken during the rising phase with C > 0.8 are consistent 
with the best fit. Those with the X-ray colour lower than 
~0.8 (and L > 0.2Le) seem to deviate from the best fitting 
line. However, the photon index and colour anti-correlation 
holds also for them, with perhaps a different slope due to 
relatively high luminosity of these hard state spectra. 

4.3 Correlation between In/ is and luminosity 

Having fitted the hard state spectra of the two sources with 
Model 2, we can now investigate the observed spectral vari- 
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Figure 4. The hard state heating-to-cooling ratio of the corona, 
£h/£s as a function of the luminosity in Eddington units for (a) 
GRO J1655-40 and (b) GX 339-4. Filled symbols - rise; open 
symbols - decay. Solid lines indicate the best fitting model to the 
data with Lbol below L cr a (see Sec. 4.3). Best fitting slopes are 
«gro = 0.60 ± 0.03 and a G x = 0.49 ± 0.15. 



ations in a physically meaningful way, by studying varia- 
tions of the hard-to-soft compactness ratio with luminosity. 
For typical coronal temperatures, the ratio of the hard-to- 
soft compactnesses determines the photon index of the hard 
X-ray continuum (e.g. Beloborodov 1999b). Consequently, 
under the hypothesis of inverse-Compton produced X-rays, 
the observed spectral shape variations should be caused by 
the lh/ls variations, which are associated mainly with the 
changes of the accretion geometry, e.g. the inner edge of the 
disc proceeding towards the black hole, or the presence of a 
coronal outflow with variable velocity. 

We present the evolution of our best fitting lh/ts values 
as a function of luminosity in Fig. [3] These plots clearly in- 
dicate that the observed decrease of the photon index, and 
the increase of the X-ray colour with increasing luminosity 
up to a few percent of the Eddington limit, is due to the 
fact that 1^ /£ a increases from a value of ~ 1 to ~ 10-15 for 
the same range of Lboi/^E, in both sources. At Eddington 
luminosity ratios higher than a few per cent the lh/ts evolu- 
tion changes. It decreases with luminosity leading to softer 
spectra (and lower values of X-ray colour). 

In order to find the critical value of luminosity at which 
the compactness ratio reaches maximum, we fit a broken 



Table 2. Coefficients in the lh/ts vs. ^bol/^E correlation. 



Source 




a 


b 




Lci-h/Le 


GRO J1655-40 


0.60 


±0.03 


328 ± 


47 


0.006 


GX 339-4 


0.49 


±0.07 


107 ± 


39 


0.01 



power law model to the lh/ts vs. I/boi/^E data of the two 
sources and summarize the fit results in Tab. 2. Similarly to 
the case of the photon index and X-ray colour correlation, 
we neglect measurements that are upper limits on lh/ts (but 
we include them in the plots), and use the average of the 
lower and upper lh/ts errorbars (from Model 2) as the error 
on lh/ts - In the case of GRO J 1655-40 we fit separately the 
data from the rise and decay and we calculate their weighted 
average, which gives Lcrit, gro = (0.0061±0.0002)Le- In the 
case of GX 339-4 we treat each outburst separately, and we 
fit the rise and decay together for outbursts ol and o3. We 
use only the decay data of outburst o2 because the data of 
the rise do not allow to put meaningful constraints on the 
critical luminosity, and a strong hysteresis effect prevents a 
joint fitting of the rise and decay data. Again, we calculate 
the weighted average of the three measurements, and we get 
L cr it,Gx = (0.011±0.002)Le. During the fits we neglect sev- 
eral data points indicated in Fig. [3]with light gray symbols. 
These observations correspond to the periods shortly be- 
fore/after the hard-to-soft/soft-to-hard spectral transition, 
when the compactness ratio, as well as the X-ray colour and 
photon index, change rapidly with luminosity. 

However, the most significant contribution to the un- 
certainty on L cr it is introduced by the uncertainty of the 
GBH distance and mass estimates. Taking into account 
these uncertainties in the case of GX 339-4 (Tab. 1) we 
derive 0.003Le < L cr it,Gx < 0.04Le, while in the case of 
GRO J1655-40 we get much tighter constrains, 0.005Le < 
L cr it,GRO < 0.007Le- Foellmi (2009) derives an upper limit 
to the GRO J 1655-40 that differs substantially from the ac- 
cepted value of 3.2 kpc. Adopting their upper limit of 2 
kpc and M = 6.3Mq we obtain Lcrit, gro = 0.002Le- How- 
ever, Caballero Garcia et al. (2007) pointed out that if GRO 
J1655-40 were indeed located so close, its companion star 
would not fill its Roche lobe. Given all the uncertainties, 
and the range of possible values for L cr it that we presented 
above, it seems quite possible that L cr it/LE is the same in 
both objects, and roughly equal to ~0.01. 

Having determined the luminosity at which the evolu- 
tion of lh/ls changes, we proceeded to determining the slope 
of the lh/ls vs. Lboi/^E positive correlation below L C rit- We 
fitted jointly the rise and decay data with luminosities be- 
low icrit.GRO and L cr it, gx, respectively, with the following 
function: lh/ls — b x (Lboi/-t/E) a - The model describes well 
the anti-correlation between lh/ls and luminosity (\ 2 /d.o.f. 
equals 1.5 and 0.7 in GRO J1655-40 and GX 339-4, respec- 
tively). The solid lines in Fig. UJi/b indicate the best fitting 
models. The best fitting values are: (Xgro = 0.60 ± 0.03, 
a G x = 0.49 ± 0.07, & G ro = 328 ± 47, fe G x = 107 ± 39 
(Tab. 2). They are consistent within their errors for the two 
sources, and their weighted means are a — 0.58 ± 0.02 and 
b = 197 ± 30. 
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Figure 5. Comparison of the hard state and quiescent state. 
Data and fits to GRO J1655-40 and GX 339-4 are as in Fig. 4 
(open/gray symbols, dashed and dotted lines, respectively). The 
filled/black symbols represent the measurements for several GBHs 
collected from the literature (see Sec. 5.5). In quiescence the lulls 
saturates at ~1.7 (solid horizontal line). 

4.4 Quiescence 

The anti-correlation between lh/ls and Eddington luminos- 
ity ratio does not extend to the quiescence. This conclu- 
sion is based on the Eddington luminosity ratios and pho- 
ton indices collected for quiescent accreting GBHs by Cor- 
bel, Tomsick & Kaaret (2006, XTE J1118+480, A0620-00, 
XTE J1550-564, GX 339-4, GRO J1655-40, V641 Sgr). We 
included also three measurements for V404 Cyg. In two 
cases, we combined the reported photon indices with the Ed- 
dington luminosity ratios calculated based on the 0.3-10 keV 
XMM-Newton unabsorbed flux (Bradley et al. 2007) and 
3-9 keV Chandra flux (Corbel, Koerding & Kaaret 2008), 
assuming M = 12M and D = 3.5 kpc (Corbel et al. 2006). 
The third quiescent observation of V404 Cyg was originally 
published by Kong et al. (2002) . However, it was claimed to 
be affected by pile-up, so we addopted the pile-up corrected 
photon index of Corbel et al. (2008), and the Eddington 
luminosity ratio as in Kong et al. (2002) . 

These quiescent measurements of the photon index are 
consistent with being constant, T q = 2.13 ± 0.06, as a func- 
tion of luminosity. In Fig. 5 we plot again the hard state data 
of GRO J1655-40 and GX 339-4 together with l h /l 3 calcu- 
lated based on the quiescent photon indices using the follow- 
ing equation, F = 2.33 x (£h/4)~ 1/6 (Beloborodov 1999b). 
We conclude that in quiescence the heating-to-cooling com- 
pactness ratio saturates at £h/£s ~ 1.7, corresponding to 

r q . 



5 DISCUSSION 

We studied two confirmed GBHs, namely GX 339-4 and 
GRO J1655-40, considering all their archival RXTE obser- 
vations until June 2007. We focused on a subset of observa- 
tions, when the photon index was smaller than 2, i.e. when 
the sources were mainly in the hard state. During these ob- 
servations, the source luminosity was lower than 10% of the 



Eddington limit in the case of GRO J1655-40, and 30% in 
GX 339-4. 

We used traditional X-ray colours and the results of 
Sobolewska et al. (2009) to determine the spectral evolution 
of each source. Our results are consistent with what has been 
observed in the past, for other GBHs and AGN as well: the 
X-ray colour increases, and the photon index hardens, with 
increasing luminosity up to the point where Lboi/LE ~ 0.01. 
Above this point, the opposite trend is observed: a decrease 
of the X-ray colour and a softening of the photon index with 
increasing luminosity. 

Within the context of thermal Comptonisation models 
for the X-ray production in accreting black holes, the photon 
index is determined by the hard-to-soft compactness ratio, 
ih/is- The main result of our work is that the observed spec- 
tral variations of the two sources, when in the hard state, can 
be explained if £h/ls and luminosity are anti-correlated for 
luminosities below L C rit ~ O.OILb, and positively correlated 
for luminosities higher than Lcrit- The anti-correlation below 
L cr it can be described by lh/ts oc L h a ol , where a = 0.58±0.02 
is the weighted mean for both objects. 

Since the compactness ratio, Ih/t-s, depends mainly on 
the geometry of the accretion flow, our results can put con- 
strains on various models that have been proposed to explain 
the hard X-ray emission in GBHs. We discuss below some 
of these constrains and implications of our results. We note 
that in the hard state spectra of GBHs the Comptonised 
emission dominates over the emission of the seed photons, 
so that Lboi oc £h- 



5.1 Hot thermal inner flow, truncated disc seeds 

The hard X-rays could be produced by Compton up- 
scattering of seed photons from a disc illuminating a hot 
inner flow. Heat conduction from thermal contact of the hot 
flow and cool disc can lead to evaporation of the inner thin 
disc, resulting in a truncated thin disc, with a radius of 
i?disc oc rh -1 / 2 (Czerny, Rozariska & Kuraszkiewicz 2004). 
We assume that the hot flow is powered by the same mass 
accretion rate that flows through the thin disc, so that its 
luminosity is given by the remaining potential energy from 
the inner edge of the truncated disc, iidisc, to the inner edge 
of the hot flow, R ln . Then (approximately) 

th ~ GMm(l/R in - l/i?disc) X 77corona/j?disc, 

where r/ corona /fjdisc is the radiative efficiency of the hot flow 
relative to that of a thin disc, and is proportional to m for 
a radiatively inefficient flow (e.g. Sharma et al. 2007), and 
roughly constant for a radiatively efficient flow. If we as- 
sume that R disc > Ri n , then l h oc mr] corona /ri disc , which 
implies that lu oc m 2 for a radiatively inefficient flow, and 
lh oc rh for a radiatively efficient flow. The truncated disc 
luminosity is L d i sc = GMrh diBC /2R <1 i BC , but only part of this 
is intercepted by the hot Comptonisation region, which we 
assume to be a sphere of radius R h . This fraction, /, varies 
with radius, and reaches maximum at the truncation radius 
of the disc, where / = (Rh/ Rdisc) 2 Since the disc lumi- 
nosity peaks at the same radius, the seed photon luminosity 
for the hot flow is 

£ s « GMm/2R disc x {R h / R disc ) 2 / tv , 
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so that 

£ h /£ s = 2tt x (R diBC /R in ) x (R diBC /R h y x n 

corona / '/disc j 

and we get that £ h /£ s oc m~ 3/2 7) corona /j7 disc oc £ h ~ 1/4 for a 
radiatively inefficient flow, and proportional to £^ Z ^ for a 
radiatively efficient flow. Either way, this predicts that £h/£s 
increases as the luminosity decreases, i.e. that the spectrum 
hardens at lower luminosities. This is due mainly to the 
rapidly decreasing solid angle subtended by the hot flow to 
the truncated disc. However, this is in marked contrast to 
the observed softening of the spectrum seen at the lowest 
luminosities, followed by the saturation in quiescence, ruling 
out this scenario as the physical mechanism for the observed 
X-ray emission at luminosities below L C rit- 

However, this could instead explain the observed be- 
haviour above L cr it, where the spectrum abruptly softens 
with increasing luminosity. Hence, there must be a transi- 
tion to another mechanism at L cr it to explain the change in 
£h/£s behaviour. 

5.2 Hot thermal inner flow, cyclo-synchrotron 

seeds 

There are also seed photons generated in the hot flow it- 
self, by cyclo-synchroton radiation of the energetic elec- 
trons. This emission is strongly self absorbed for ther- 
mal electrons, so it peaks at the self absorption frequency 
zAjsa, given where the cyclo-synchrotron emissivity equals 
that of a blackbody. The full equation for this is complex 
(Narayan & Yi 1995). However, Wardziriski & Zdziarski 
(2000) give a numerical approximation for this frequency, 
i/ ssa tx 0°- 95 r o ' O5 f? ' 91 , where and r are the dimension- 
less temperature and optical depth of the electrons in the 
hot flow, and B is the magnetic field. At the observed tem- 
peratures, z/ ssa should be on the Rayleigh- Jeans tail of the 
blackbody, so the flux at this point is oc v 2 s ^9. Thus the 
seed photon luminosity l s oc z/g Sa #. The hard luminosity 
is as before, with £h = GMm/2Ri n x 77corona/?7disc, hence 
£ h /£ s oc V(0 a - 85 r°- 15 B 2 - 73 ). 

We get a further constraint as £yj£ s is also set by 9 and 
r. Pietrini & Krolik (f995) give an approximation for this of 
O.f (£ h /£s) 1/4 = Br. This gives l h /l s oc £^ 51 t 1a9 (r / B 2 ) - 70 . 
In a radiatively inefficient flow, the gas pressure is set by 
the ion temperature, which remains approximately constant 
around the virial temperature. Hence t/B 2 is constant if the 
magnetic pressure is a constant fraction of the gas pressure. 
Hence £ h /£ s oc < k ' 51 r us oc £ h 1A if r oc m oc £^' 2 . This 
is somewhat faster than the observed £h/£s oc £®"'~ ' 6 , but 
given the level of approximation this at least describes the 
observed trend below L cr it for the spectrum to soften as the 
luminosity decreases. 

However, this alone does not give a framework for the 
change in behaviour at higher luminosities, where the spec- 
trum dramatically softens with increasing luminosity. Thus 
neither seed photons from the disc nor self produced syn- 
chrotron self Compton photons from the flow can explain 
the full range of observed behaviour, though the combina- 
tion of the two (seed photons from the disc changing to seed 
photons from cyclo-synchrotron as luminosity decreases) can 
explain the observations. 

The saturation observed in quiescence is not expected 



if the electron distribution in the hot flow is thermal. How- 
ever, it is likely that the electron distribution is not a pure 
Maxwellian. Evidence for the presence of a non-thermal high 
energy tail in the electron distribution was found in the hard 
state of several black hole binaries (McConnell et al. 2002; 
Wardzihski et al. 2002; Cadolle Bel et al. 2006; Joinet et 
al. 2007; Droulans et al. 2010). Consequently, in quiescence, 
while the plasma becomes optically thin (due to decreas- 
ing L/Le), the Comptonised component becomes dominated 
by the contribution arising through scattering on the non- 
thermal high-energy electrons, and the evolution of T (and 
hence th/£s) depends on the microphysics of particle acceler- 
ation in the corona. Particle acceleration with a slope that 
does not depend on L/Le would result in a saturation of 
£ h /£ s . 

5.3 Outflowing hot thermal corona, untruncated 
disc 

An outflowing hot corona above an untruncated disc is an 
alternative geometry to explain the observed hard spectra 
seen in the hard state (Beloborodov 1999a; Malzac et al. 
2001). Here the seed photons for Compton scattering are 
from the disc, but are suppressed by relativistic beaming 
from the mildly outflowing corona. 

In the original model of Beloborodov (1999a), the in- 
trinsic disc emission is assumed to be negligible compared 
to reprocessed radiation produced by the disc illumination 
by the X-ray source. Then Ih/ls depends mainly on the out- 
flow velocity (as well as a few geometric parameters) . As the 
velocity is increases, beaming reduces the illumination of the 
disc and the soft photon flux returning to the outflow. There- 
fore, as the velocity increases, the ratio £h/£s increases, and 
the reflected fraction decreases. We also note that in the 
outflowing corona model, the luminosity of the corona is 
amplified by relativistic beaming and this affects the shape 
of the observed correlations between £h /£ s and luminosity. 

The observed spectral softening (and correlated increase 
in reflected fraction: e.g. Gilfanov, Churazov & Revnivt- 
sev 1999; Ibragimov et al. 2005) seen above Z/ cr it can be 
explained in this model if the outflow velocity (and hence 
th/£s) decreases as luminosity increases. The change in be- 
haviour below -Lcrit then requires that the outflow velocity 
decreases as luminosity decreases. This makes a clear pre- 
diction that at the lowest luminosities where the beaming 
is negligible, the continuum should be accompanied by a 
substantial reflected emission from an untruncated disc and 
£h/£s should saturate, as observed. 

Alternatively, the softening of the spectrum below i C rit 
could also be caused by an increase in intrinsic emission 
from the disc such that this dominates over the reprocessed 
emission as the source of seed photons. The outflow veloc- 
ity could then remain high, or even continue to increase as 
luminosity decreases. This predicts that the reflected frac- 
tion remains small, but that the observed fraction of power 
dissipated in the disc should increase. 

5.4 Non-thermal jet, untruncated disc 

Another type of emission suggested to explain the hard X- 
ray spectrum is direct non-thermal synchrotron radiation 
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from a jet (Markoff, Nowak & Wilms 2005). Both, thermal 
Comptonisation and jet models can reproduce the shape and 
observed energy range of the spectral cut-offs seen in the 
hard states (e.g. Zdziarski et al. 2003; Markoff & Nowak 
2004). In several black hole binaries the cut-off energy was 
reported to anti-correlate with luminosity in a bright hard 
state (e.g. Wardziriski et al. 2002; Rodriguez, Corbel & Tom- 
sick 2003; Miyakawa et al. 2008; Joinet, Kalemci & Sen- 
ziani 2008; Motta, Belloni & Homan 2009). This can be ex- 
plained rather naturally in terms of thermal Comptonisation 
models by enhanced cooling by the disc photons. Thus, it 
seems likely that thermal Comptonisation does dominate at 
Lboi/Ls > 0.01. It is then possible that the transition in 
spectral behaviour seen at L cr it marks the point at which 
non-thermal emission from the jet starts to dominate over 
thermal Comptonisation from the hot flow (Russell et al. 
2010), and theoretical models of jet dominated accretion 
flows (e.g. Falcke & Markoff 2000; Markoff et al. 2001; Yuan, 
Markoff & Falcke 2002) certainly predict that the kinetic 
power of the jet exceeds the radiative luminosity (see also 
Malzac, Merloni & Fabian 2004; Malzac, Belmont & Fabian 
2009). However, the photon index of the X-ray emission 
then depends on the photon index of the non-thermal elec- 
trons within the jet, so to explain the spectral softening 
with decreasing luminosity requires that the electron heat- 
ing/cooling processes result in a steeper electron distribu- 
tion. There are currently no real constraints on this, due to 
the lack of knowledge of jet physics. 

In quiescence, similarly to the case of non-thermal 
Comptonisation dominating the X-ray in the scenario in- 
volving thermal inner flow and cyclo-synchrotron seeds, the 
saturation of £h/£s can be explained in the jet model if the 
slope of accelerated non-thermal particles in the jet is inde- 
pendent on L/Le- 



6 CONCLUSIONS 

We show that in the hard state of GBHs the compactness 
ratio £h/is reaches maximum at L cr it ~ O.OILe- For Lboi < 
icrit, we find that £h/£s = bx(L\ }0 \/ Le) a , where a = 0.58± 
0.02 and b = 197 ± 30 are the weighted slope and intercept 
based on the fits to the data of GRO J1655-40 and GX 339- 
4. We suggest that the change of the behaviour in the £h/£s 
vs. luminosity relation is caused by a change in the X-ray 
emission mechanism for the hard state spectra of GBHs. The 
relation does not extend to the quiescent state. Instead, a 
saturation of £ h /£ s is observed at a value of ~ 1.7. 

The observed hard state evolution of £u /£s can be ex- 
plained if the seed photons are cyclo-synchrotron photons 
up to Lcrit, above which they are replaced with the trun- 
cated disc photons. Alternatively, the observed variability of 
£h/£s can be explained in the case of an outflowing corona 
if the outflow velocity increases with luminosity up to Lcrit 
and then decreases with increasing luminosity, or if the seed 
intrinsic disc photons are replaced with the reprocessed pho- 
tons above Lcrit- Finally, in the scenario involving an X-ray 
emitting jet, the observed transition in £h/£s behaviour may 
mark the point at which the non-thermal jet emission starts 
dominating over the thermal Comptonisation emission. 

Observed saturation of £h/£s in quiescence can be ex- 
plained in terms of an untruncated disc model with cyclo- 



synchrotron seed (and in the jet model) if the non-thermal 
particles in the corona (jet) are accelerated with a slope that 
is independent on L/Le- The outflowing corona model natu- 
rally predicts the quiescent saturation of £h /£s if the outflow 
velocity decreases with luminosity below L cr ; t . 
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